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Abstract 

We prove the renormalizability of the Curci-Ferrari model with and without 
auxiliary fields using BRST methods. In both cases we find 5 Z factors 
instead of 3. We verify our results by explicit one loop calculations. We 
determine a set of generators for the "physical states" , many of which have 
negative norm. Supersymmetrization is considered. 
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Recently, there has been renewed interest in consistent actions for massive 
vector bosons without Higgs fields. In particular, Periwal has reanalyzed 
the action of Curci and Ferrari (0, see also and ||) for massive vector 
bosons and claimed that it is both renormalizable and unitary If true, this 
would be astonishing since a great deal of work over the past 30 years has left 
only the Higgs mechanism as a means of giving vector bosons a mass while 
preserving renormalizability and unitarity. Another recent suggestion is to 
use topological field theories]^]. In fact, in his pioneering work in the 60's, 
Veltman|| started with massive Yang-Mills theories coupled to free scalars, 
but after a field redefinition which made the scalars seem to interact, the re- 
quirement of renormalizability forced him to drop certain terms, thus ending 
up with the Higgs model. Other studies]^]] which required that tree graphs do 
not grow too fast with energy in order to obtain one-loop renormalizability 
confirmed these results in the 70's. 

A consistent model for massive vector bosons would not yet be an alter- 
native to Higgs fields since it is not immediately clear how to give a gauge- 
invariant mass to the fermions of the standard model. It would, however, 
allow infrared regularization of QCD and of supersymmetric models (in par- 
ticular in superspace) both of which are plagued by serious infrared problems. 

The model of Curci and Ferrari (CF) has both a mass term for the vector 
bosons and a gauge fixing term. As a result the propagator has k^k v /k 2 
terms instead of a non-renormalizable k^kv/m? term, and is therefore power- 
counting renormalizable. It also has a BRST (and an anti-BRST) symmetry, 
but here the analogy with Higgs models stops: the BRST operator is neither 
nilpotent, nor can it be made nilpotent by introducing a BRST auxiliary 
field. Despite much work in the past, it is not clear whether this model can 
be obtained as a suitable limit of a Higgs model. As a result of the non- 
nilpotency, the usual "rr" Ward identity of non-abelian gauge theories || is 
modified by an extra term, and the issue of renormalizability requires deeper 
study. In this letter we first study the renormalizability of the CF model, 
both with auxiliary field and without, and then we shall come back to the 
issue of unitarity. We find that the theory is renormalizable, as already 
found by Curci and Ferrari using different, more cumbersome, methods, and 
by Periwal, but we find 5 Z factors instead of 3, as claimed by Periwal. Then 
we determine the physical states, extending Ojima's workH. Many of these 
states have, for arbitrary values of the parameters of the theory, a negative 
norm, and from this we conclude that the model is not unitary. Finally we 
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briefly discuss supersymmetrization of the model. 
The CF action is given by 

S = Sym + S m + 5gf, (1) 

where Sym is the Yang-Mills action — \(F® U ) 2 , S m contains a mass for the 
vector bosons and the ghosts c a and antighosts b a , while contains terms 
which resemble the Faddeev-Popov ghost action and the gauge fixing term 

S m = -^m 2 ((^) 2 + 2^ a c a ), (2) 

Sgf = ~b a (d»D^ + D^)c a -±(d»Al) 2 + ±g 2 Z(bxc) 2 , (3) 

where b x c = f a bc b b c c . Both S m and S g f are separately BRST invariant; A ^ 
and c a transform as usual (5A^ = D^cA and 5c = |gcxcA; we omit the BRST 
parameter A) while the BRST law of b a can be found by requiring invariance 
of S m . It reads 5b = — 1<9- A+ \gb x c. Since 55b is no n- vanishing, the action 
Sgf can be found by assuming that for vanishing mass 55b is proportional to 
the b field equation, and integrating the latter. 

We first perform the analysis of renormalizability with a BRST auxiliary 
field A a present. If we define 5b = A but 5X = —m?c (rather than 5\ = as 
nilpotency would require) the action S — Sb is BRST invariant, where 

S = Sym + S m + + Sb, (4) 
S^ = S gf + ^(A + ^-^xc)) 2 (5) 

= b a d»D,c a + ^{b x c) 2 - A;d,X a + l -i\ 2 - ^gi\ ■ (b x c) 
= 5[b a (^X a + dA a - \giib x c) a )} + \m 2 ib a c\ 

S B = J[K ■ DC + L-gcX c + M ■ A + N{-m 2 c)]. (6) 

Following Zinn- Justin and B. Lee||, we couple the BRST variations to ex- 
ternal sources. We add source terms for the fields, S s = J[JA + lc + mb + nX], 
and obtain then the Ward identity for the effective action T = W — S s 

d d d d d d d d 
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As a check one may verify that T = S satisfies this equation. In the renor- 
malized theory, all Z factors in this relation should amount only to an overall 
rescaling. Hence, defining = Z^ 2 A^(ren) etc., Kg = Z^ 2 Kg (Ten) etc., 
£ = Z^(ren), m 2 = Z m 2m 2 (ren), g = Z g g(ren), we assume inductively the 
relations 

^1/2^1/2 _ 7I/27I/2 _ 7 l/2 7 l/2 _ „l/2„l/2 _ 7 -l 7 -l/2 7 l/2 _ 7 -l/2 7 l/2 
Z-3 Zj k — Zj c Zj l — Zj h Zj m — Zj x Zj n — A m 2^ N 6 L — 6 M 6 N ■ 

(8) 

Since the action has vanishing ghost number, we expect that only the product 
(Z b Z c fl 2 = Z z will be fixed. 

Since V contains only one-particle irreducible graphs, it is independent 
of M and N, so we drop the term — M-^T term. Assuming (l — 1) loop 
finiteness and the scaling hypothesis in the /-loop divergences T(div) 
must satisfy Qr(div) = 0, where 

with x 1 = {A, — L, — M, A} and &i = {K,c,b, N}. We decompose Q and 
T(div) into terms without m 2 and terms proportional to m 2 

q = g(o) + m 2gd). r(div) = r(°) + m 2 r( 1 ) (10) 

where S{m 2 = 0) = SVm + 5rf + K,L,M terms. Again we drop the term 



ic-J^- r We must then solve 

Q(°)r<°> = 0; g (0) r« + Q«r(°) = o ; g«r« = o. (13) 

Since Q 1 - -* is nilpotent (the term N(—m 2 c) is no longer present in S(m 2 = 0)) 
we have 

r (0) =(3S YM + Q {0) [X + Y], (14) 

where X and Y denote the most general terms which are M, N independent 
and M and/or N dependent, respectively. Clearly, 

X = a\K ■ A + a 2 6 • dA + a 3 L • c + a^gb ■ b x c + a^b ■ A, (15) 
y = a 6 gN ■ b x c + c^iV • X x c + a 8 iV • A + a 9 M ■ b 

+a 10 M • N + a n N ■ dA. (16) 
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Q^X is already M and N independent, but requiring Q^Y to be M, N 
independent leaves only 

Y = a 6 (M -b + N ■ X). (17) 
The most general form of is 

r« = |( 7l A 2 + 72 ^-c), (18) 

which already satisfies Q^T^ = 0. Hence we only need to solve Q^F^ + 
q(!);T(o) — o. This yields the following three relations in the nine coefficients 
. . . ,a 6 ,7i,7 2 

2«i — a 2 — a 3 + a 6 + 2 7i = (19) 

2a A + g£a 6 + ^ 72 = (20) 

2a 5 - 2^6 - £72 = 0. (21) 

However, in T(div) only 5 combinations of parameters occur because the 
combination Q(°\b ■ N) in X + Y is obviously annihilated by Q(°\ 
The final form of the divergences is 

r(div) = ps YM + 

+ a^A^-Syu - m 2 A 2 + \ ■ dA - gdb ■ A x c - K ■ dc] 

+ a[-^m 2 A 2 + X-dA + bd 2 c-gdb-Ax c] 

+ l2 [m 2 &c - - A ig 2 {b x c) 2 - ^A 2 + \tg\ • b x c] 

+ a^m 2 A 2 + \ig 2 {b x c) 2 + bd 2 c -gdb-Axc 

-i^A • 6 x c + K ■ dc + L^yc x c + K ■ gAx c\. (22) 

As a check we have verified that each of these terms is annihilated by Q. In 
the process one needs the identities 

(^S YU ) ■ dc = 5(A-^S YM ); (b x cf = ~(b x b) ■ (c x c); (23) 

(ylxc)xc=^x(cxc)i c-bx c = b ■ cx c. 
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Expanding the Z factors as Z g = 1 — z g , Z 3 = 1 — za etc., we find 



z g = --(3; z A = 2«i + P; z m 2 = -(3 - a 3 + a; z i = [3 - 72 - 2a 

z x = 2a- (3; ^{z b + z c ) = a + a 3 ; (z K + z c ) = 2a 3 - 2a x \ 

X -z L + z c = a 3 + ~/3. (24) 

As expected from ghost number conservation, only Zb+z c , zk + z c and z^+2z c 
occur. Since M and iV do not contribute to V except at the tree level, we 
further get zm = —z\ and z N = —2z m 2 — z c . One can easily check that 
these Z factors satisfy the scaling hypothesis (H). This completes the proof 
by induction of the renormalizability of the model with the auxiliary field 
present. There are clearly 5 Z factors. 

The action has an infinitesimal U(l) invariance 5c = b,5b = — c, 5X = 
\gb x b — \gc x c, leading to a finite symmetry 

c ^ b,b — > — c, A — >• A — gb x c. (25) 

Under this symmetry C m , £ g f and (A + jd ■ A — \g{b x c)) are separately 
invariant. However, T(div) does not have this symmetry even at K = L = 0, 
nor should it, since 5 A is nonlinear in fields and hence is modified at the 
quantum level. Since S + T(div) can be written as the renormalized S, it is 
clear that the renormalized transformation rule for 5X keeps S' + r(div) U(l) 
invariant. Requiring (erroneously) that only T(div) be U(l) invariant would 
lead to a 3 = 0, but we keep a 3 . 

Eliminating A by substituting the A field equation A = — 1<9 • A + \g(b x 
c) — |M in S and in T(div) we find 

S = Sym + S m + S'gf 

+ j[K-DC + L l -gc x c + M{--d ■ A + l -g(b x c)) - j-M 2 ], (26) 



and 



r(div, no A) = (3Sym + 



r\ 11 1 

+ ai [A—S Y M - m 2 A 2 - -{dAf + -gA^f - K ■ dc - -M ■ dA] 
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+ a[-\m 2 A 2 - -AdAf + bd 2 c + \gA»f - -M ■ OA] 
2 4 24 

+ 72 [m 2 e&c - kg 2 (b x c) 2 - hdAf --M-dA- ±-M 2 } 
° 24 4 24 

+ a 3 [^m 2 A 2 + i^ 2 (6 x c) 2 + bd 2 c + \gA.f 

1 1 
+K ■ dc + L-gc x c + K ■ gA x c+ -gM ■ b x c], (27) 

where = d^b x c — 6 x <9 M c. Notice that A only couples to the ghosts 
through the U(l) current. All the terms in r(div, no A) are U(l) invariant 
when K = L = M = 0. One can easily check that the same Z-factors still 
render the model finite. 

Furthermore, as an additional check of our results, we began with the 
action (p6|) without auxiliary fields and performed an analysis similar to that 
given above. This is possible because the term — (1/24')M 2 plays the role of 
the usual subtraction S — » S — «Sg x that one makes when solving the BRST 
cohomology in the usual case of non-abelian gauge theories 0. Of course, the 
results are exactly as given in fl2"7p. 

Periwal in has shown that the same model is renormalizable, however 
with two fewer parameters as we now discuss. Comparing our results to those 
of Periwal, we find that his results are a subcase of ours in the following sense 

= Z 3 3 Z 3 l ^ 2 Z g 3 ; Z m 2 = Z^Z 3 Zg. (28) 

The conditions in (p8[) imply 

72 = 03; a + ai + 2a 3 = (29) 

Then he finds 

%k = Z 3 3 Z 3 2 Z g 4 ; Z L = Z M = Z 3 4 Z 3 l Z g 4 (30) 

which agrees with our results if we use (|29"D. To check that the relations in 
(|29f) are not due to a symmetry which we overlooked, we made an analysis 
of divergences at the one loop level. 

Table 1 gives the result of one loop calculations for the divergences indi- 
cated, followed by combinations of parameters which get fixed in this way. 
The number C 2 is defined by f ap q f b q p = —C 2 S^, and D denotes the standard 
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g*C 2 DJbdic(-l + t) 


a + 0:3 


g 2 C 2 DJb^m 2 c(^) 


72 


g 2 C 2 DJKdc(f) 




g 2 C 2 DJgK-Axc(i) 


«3 


g 2 C 2 Dj\gL-cxc{^) 




g 2 C 2 D j\gA- (db x c-bx dc) (§) 


a + «i + «3 


g 2 C 2 DJm 2 A 2 (-§ - |) 


"«i - f + f 


^ 2J D/(^) 2 (i|-|) 


-f-ai 


g 2 C 2 DJ(d-A) 2 (^ + i) 


§ + «i-?- 


g 2 C 2 D J =±M ■ OA (0) 


a + «i + 72 


g 2 C 2 Dj^M 2 (^) 


72 


g 2 C 2 Dj\gM-bxc(^\ 


«3 



(31) 



Table 1: One loop divergences. The column on the right indicates which 
parameters are fixed by the corresponding divergence. 
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\g 2 iC 2 D 
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g*£C 2 D 



Figure 1. 



divergence / d 4 k{2^) A (k 2 + m 2 ) 2 . From these results we conclude that at 
the one loop level the divergences are given, up to an overall factor g 2 C 2 D 
by ^ ^ ^ ^ 

/5 = -y, ai = -(3 + 0, « = 72 = a 3 = -f (32) 

Clearly, the conditions in fl2~9|) disagree with the values in (|32|). In particular 
the quantity gA/m 2 ^ does renormalize. 

The simplest check that the relations (|29|) do not hold, is given by the 
graphs in figure 1. The first one yields 72, while the sum of the last two 
graphs yield 0:3. Clearly, at the one loop level 

"3 = -72, (33) 

which does not agree with (29). 

Although the model is renormalizable, it does not seem to be unitary. We 
perform a Hamiltonian analysis to determine the spectrum of the theory. We 
assume the usual relation between Heisenberg fields and in- and out- states 0. 
We eliminate and p{A^) from the Hamiltonian by using the second class 
constraints p(X) — A = and p(A Q ) = 0. A canonical transformation with 
generator 

F 2 {q, P)= J d 3 x{(A? + -Ld^PiiU) + A a P a (A)] (34) 

decouples the A-terms from the A- terms in the Hamiltonian. The Hamilto- 
nian density in terms of the new asymptotic fields Uf and A a reads 

H = \p\8 ij -^)pl-\u i {d 2 5 ij -d i d j -m% j )U j 
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~m 2 P 2 + 7^A(<9 2 - m 2 OA - p(c)p(6) - &(9 2 - m 2 £)c, (35) 



where <9 2 = d l di. Note the unphysical sign in front of the the A terms. All 
field equations are of Klein-Gordon type with m 2 for [/, and Pi, and £m 2 for 
A and P(A). 

The corresponding modes satisfy the commutation relations 

[<(£), ajflt] = 5 ab (5 tJ + ^)5(k - f), (36) 

[<(£), a^^] = -5 afe 5(£ - T). (37) 

For the ghost/antighost system one finds the usual {c a (k), b b {ly} = 5 ab 5{k — 
I). The BRST tranformations read 

5di(k) = 0; 5b{k) = ma A {k); 5a\(k) = mc(k); 8c(k) = 0. (38) 

These results are the canonical counterpart of Ojima's analysis in 0, where 
he shifted into U" = A a + <9 M A a /m 2 , after which and A a decouple in 
configuration space. 

Defining "physical states" to be BRST and anti-BRST invariant, all ai(k)^ 
generate physical states. However, in contrast to the usual case, there are 
further physical states made up from ghosts, antighosts and A oscillators, 
some of which have negative (and ^-independent!) norm. One might expect 
this since Q is not nilpotent, as a result of which the unphysical states do not 
form Kugo-Ojima quartets. We restrict our analysis here to the case with 
the auxiliary field present, but similar remarks apply also to the formulation 
without the auxiliary field. In that case there is a particularly simple ghost- 
dependent observable, namely the mass-term in the action. 

To find the states in the formulation with the auxiliary field, we note that 
the anti-BRST transformations read (from now on we take m — 1) 

5ai(k) = 0; 5c(k) = -a A {k); 5a A (k) = b(k); 5b(k) = 0. (39) 

Clearly, 5,5 and {5,5} form a graded version of sl(2), and (b, A, c) form a 
spin-1 representation. The BRST and anti-BRST invariant states are the 
singlets of this graded si (2). To find the physical states we must determine 
the singlets in tensor products of the spin-1 representations. This problem 
has been solved for the usual (ungraded) si (2) case by Weyl[|I0|j. 
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All invariants are built either from inner products (vi,V2), or from deter- 
minants det(vi, v 2 , V3), with Vi three-component spin-one vectors. Adapted to 
our case, this yields the following (overcomplete) set of generators of physical 
states 

(40) 
(41) 



(42) 
(43) 

and q = c ai (ki)^. To obtain a physical state, we still need to contract a 
product of these generators with a suitable tensor, so as to make the final 
result invariant under global gauge transformations. Under the latter, A r (ki) 
transforms as (T ai ) r+2 , A ri > r *{k x , k 2 ) as (T ai ) ri+1 (T a2 ) r2+1 , etc. None of the 
states made up out of these generators are of the form 5X + SY. They 
have vanishing ghost number and are invariant under the linearized U(l) 
symmetry b — > c, c — > —b. The state A°'°(k, t)\0) was found by Ojima J4|, 
who observed it has a negative norm. In fact, many of the physical states in 
this sector of the Hilbert space have negative, ^-independent norm, as can 
be seen from the examples 

K(fc)|o>l 2 = 

|A n ' r2 (A;i, A; 2 )|0)| 2 = 

We can conclude that the CF model is renormalizable but not unitary for 
any value of £. Of course, this argument would break down if the relation 
between Heisenberg fields and in- and out- states || is no longer valid in this 
model. This might for instance happen if the theory only makes sense if it is 
strongly coupled, or has bound states, but then it is not clear to what extent 
a perturbative analysis can be trusted. 
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A r (h) = (r + 2)\ r 1 b 1 c 1 -X[ +2 
A r ^{k l7 k 2 ) = H X H 2 - Ai 1 A2 2 (&ic 2 + b 2 c 1 ) 
A r ^(h,k 2 ,k 3 ) = H x \^\ r i{b 2 ^ + hc 2 ) 

+H 2 X^XY(b 3 c 1 + b 1 c 3 ) 
+H^\ r i(b l c 2 + b 2 c 1 ) 
—2HiH 2 H 3 

where 

Hi = A[ l+1 - riXl l '\c u 
the r, are non-negative integers, k\^k 2 ^ k 3 , and Aj = a e ^(ki)^, b^ 



(-l) r+1 (r + 2)r! (44) 
(-l) ri+r2+1 r 1 !r 2 ! (45) 
2(-l) ri+r2+r3 r 1 !r 2 !r 3 ! (46) 



Although the model is not unitary, it might be useful as a regulariza- 
tion scheme for infrared divergences. In particular, in superspace, where 
dimensional regularization is incompatible with supersymmetry, this scheme 
may finally resolve long-standing problems concerning infrared divergences. 
To supersymmetrize the model, one might start from the observation that 
the double BRST variation of the antighost is proportional to the antighost 
field equation for vanishing mass. Alternatively one might seek an action 
of the form S gf = Q (0) J d 8 zTrb[X + F] + h.c. and try to solve Q {l) S gf + 
Q^°'S m = for F. In the former case one finds the BRST transforma- 
tion rule for the chiral antighost b from the invariance of a mass term 
S m = (m 2 /g 2 )Jd 8 zTre v + (m/g 2 )(J d 6 zTrbc + h.c.) using the BRST trans- 
formations 5e = —ce v + e v c, 5c = — c 2 (here 5 is an antiderivation). The 
result is 5b = -mD 2 e v - (l/2){6,c}. One can then write a BRST invariant 
action 



S — SsYM + Sm + 

1 r , m 



1 

7i 



J d 8 z{-b{-ce v + e v c) - ^{b, c}e v ) 

(47) 



+ d 6 z(-^{b,c} 2 -^(D 2 e v ) 2 ) + h.c 
J 8m 2 

We note that this action does not have a gauge fixing of the form (d ■ A) 2 
in the bosonic sector. Further work on the supersymmetrization of the CF 
model and its applications is in progress [ITT]. 
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